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ABSTRACT 

The evolution of symmetric extragalactic radio sources can be characterized 
by four distinct growth stages of the radio luminosity versus source size of the 
source. The interaction of the jet with the ambient medium results in the forma- 
tion and evolution of sources with non-standard (flaring) morphology In addi- 
tion, cessation or restarting of the jet power and obstruction of the jet will also 
result in distinct morphological structures. The radio source population may 
thus be classified in morphological types that indicate the prevailing physical 
processes. Compact Symmetric Objects (CSOs) occupy the earliest evolutionary 
phase of symmetric radio sources and their dynamical behavior is fundamental 
for any further evolution. Analysis of CSO dynamics is presented for a sample of 
24 CSOs with known redshift and hotspot separation velocity and with a large 
range of radio power. Observables such as radio power, separation between two 
hotspots, hotspot separation velocity and kinematic age of the source, are found 
to be generally consistent with the self-similar predictions for individual sources 
that reflect the varying density structure of the ambient ISM. Individual sources 
behave different from the group as a whole. The age and size statistics confirm 
that a large fraction of CSOs does not evolve into extended doubles. 

Subject headings: galaxies: active - galaxies: jets - galaxies: evolution 
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1. Introduction 



The morphological shape of all ex- 
tragalactic radio sources is characterized 
by distinct structural components that 
may be more or less prominent. The 
jet-hotspot-lobe-cocoon structure of the 
largest classical double sources has first 
been classified in two distinct populations 
dFanaroff fc Riley 1974]) . FRI sources 
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have lower radio powers and less well- 
confined outer lobes energize by a lossy 
jet and hotspots that are at less than 50 
percent of the total extent of the source. 
Their outer structure becomes a flaring, 
meandering plume that is shaped by the 
motion of the galaxy or by internal fluid 
instabilities. A typical FRI source 3C 
31 has a relatively prominent inner jet 
and filamentary plumes extending 300 kpc 
dLaing et al. 20080 . In contrast, FRII 
sources have higher radio luminosity and 
well-defined outer lobes that surround 
faint jets and very prominent hotspots at 
more than 50 percent of the source ex- 
tent. FRII sources also exhibit extended 
cocoon structures in low-frequency inter- 
ferometric images. An archetypal FRII 
source, Cygnus A exhibits two symmet- 
ric edge-brightened lobes spanning over 
160 kpc, a prominent hotspot at the lead- 
ing edge of each lobe, but with a weak 
flat spectrum core and intervening jets 
fICarilli fc Barthel 19961) . This FR clas- 
sification reveals a morphological sequence 
that relates to the ability of a jet to trans- 
port and deposit momentum and energy 
at the leading edge of the lobe. 

These large-scale Symmetric Objects 
(LSO) form an extended family together 
with Medium-sized Symmetric Objects 
(MSO) and the smallest Compact Sym- 
metric Objects (CSOs), which represent 
the earliest development phase for all dou- 
ble sources. Scaling up in age and size 
from CSOs, the MSO group comprises of 
the GHz-Peaked Spectrum (GPS) sources 
fIFanti et al. 19901 IQ'Dea 19911) (review 
IQ'Dea 1998| IFanti 20090 . the Compact 
Steep-Spectrum (CSS) sources (IPeacock &: Wall 
(review IFanti fc Fanti 19941 IFanti 20090 . 



A further morphological group of core-jet 
Blazar sources with jets pointing at the ob- 
server includes BLLac Objects, High Po- 
larization Quasars and Optically Violent 
Variable Quasars. 

The structural and spectral character- 
istics of radio sources are determined by 
the power of the source, the local environ- 
ment of the host galaxy, and the evolu- 
tionary age. In this respect, the compact 
and supposedly young CSOs (and GPSs) 
are critical for understanding the rest of 
the radio source population. Their com- 
pactness has been attributed to two dis- 
tinct scenarios: (1) the youth scenario (e.g. 
Phillips fc Mutel 1982[ IFanti et al. 19951) 



suggests that they are in an early evo- 
lutionary stage, and may continue to 
grow to Mpc-scale extended radio sources, 
and (2) the frustration scenario (e.g. 
van Breugel et al. 1984[IO'Dea 199HICarvalho 19941 



ICarvalho 19980 suggests that their growth 
is retarded (stagnated) by the dense Inter- 
stellar Medium (ISM) within the host 
galaxy. Both scenarios may apply for 
CSOs (and GPSs) as they relate directly 
to the power level and the duration of 
the nuclear activity. Young sources with 
persistent long-term nuclear activity will 
continue to grow and eventually become 
LSO doubles. Other young sources may 
stagnate because of intermittent behav- 
ior of the nuclear activity and the jet 
power at any time during their evolution 
( Reynolds fc Begelman 1997 Kunert-Bajraszewska et al. 2* 

The complicated fluid and radiative be- 
havior of extra-galactic sources driven by 
the jet outflows has been modeled using 
self- similarity of the flow pattern as it 
1 Selves and moves away from the host 



galaxy into regions with lower ambient 
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densities. Using such simplified mod- 
els for well-defined doubles, the radio lu- 
minosity and spectral characteristics as 
a function of density gradient and the 
size of the source have been considered 
by a number of authors (ICarvalho 1994[ 
IKaiser fc Alexander 19971 ICarvalho 19981 
IKaiser fc Best 2007j) . The advancing mo- 
tion of the lobes and hotspots is fully de- 
termined by the momentum carried by the 
jets, the locations of shocks in the jets, 
hotspots and lobes, and the density struc- 
tures of the ISM of the host galaxy and 
the Inter-galactic Medium (IGM). On the 
other hand, the morphological radio struc- 
ture of FRI-like sources is very different 
from the structure of FRII-like sources. 
After the formation of a weak hotspot at 
a certain distance from the nuclear engine, 
the flow is conically confined and forms 
a flaring (and meandering) structure ex- 
tending for a large distance, as modeled 
for the well-known FRI source 3C31 (see 
Wang et al.~2009| ). 

The current investigation considers the 
observable properties of a well-defined 
sample of CSO sources, in order to deter- 
mine their structural evolution with time. 
The observable properties of CSO are im- 
portant to understand the early evolution 
of radio sources and how they evolve into 
the larger population of extra-galactic ra- 
dio sources. Earlier modeling using self- 
similarity of the shape of the FRII-like 
radio sources assumed a constant separa- 
tion velocity of the radio hotspot. How- 
ever, analytic modeling of the expansion of 
hotspot and cocoon indicates that the ra- 
dio sources experience a deceleration from 
CSO to MSO phases, and an acceleration 
process from MSO to LSO-FRII phases 



flKino fc Kawakatu 20051 IKawakatu fc Kino 2006]) . 
While such a trend would also follow from 
simple dynamical analysis and from the 
analogies with fluid dynamics, the hotspot 
separation velocity will be one of the ob- 
servables in the study the earliest CSO 
stage. 

The radio morphology of all mem- 
bers of the family of radio sources repre- 
sents the evolutionary stage of the source 
but also reveal symptoms of survival or 
demise. In order to connect the early evo- 
lution of CSOs with the later evolution of 
the radio source family, the physical pro- 
cesses that determine the morphology and 
the radio properties will be used to fur- 
ther classify the evolutionary phenomenol- 
ogy. Seven morphological types within the 
CSO, MSO, LSO family are used to in- 
corporate and describe what may happen 
to extra-galactic radio sources during their 
life time. 

This paper presents the following com- 
ponents: a description of the dynamical 
evolution of extra-galactic radio sources 
and the place of CSO sources within a 
morphological classification scheme (Sec- 
tion 2), a description of the observational 
data and the CSO sample used for this 
work (Section 3), the observed dynamic 
changes of CSO sources (Section 4), the 
statistics of the CSO sample (Section 5), 
and conclusions (Section 6). 

2. Evolution of the Radio Structure 

The variation of the radio power P rac i 
versus the total extent of the source D 
for the whole family of extra-galactic radio 
sources serves to further understand their 
evolutionary characteristics (Fig. [I]). Such 
diagrams have been presented earlier in 
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the literature (e.g. IQ'Dea fc Baum 1997| 
Kunert-Bajraszewska et al. 2010|). The 



current (P ra d ~ D) diagram contains simi- 
lar but also different groups of data points 
than earlier versions, particularly for the 
region of early development of CSO ra- 
dio sources. A sample of low-power GPS 
sources (Ide Vries et al. 20090 has been in- 
cluded in the bottom-left corner of the 
diagram, as likely precursors of the low- 
power MSO sources (? ). In addi- 
tion, there are more data points for the 
group of low-power FRIs (< 10 24 W/Hz) 
dLaing et al. 1983[IFanti et al. 19871 IFanti et al. 



diative properties of the hotspots and the 
well-confined lobes may thus be derived 
from modeling (IKaiser fc Alexander 19971 
IKawakatu fc Kino 2006[IKaiser fc Best 20071) . 
although the assumption of self-similarity 
limits the interpretation of some observ- 
able parameters. 

The density gradient along the path of 
the jet (i.e., in the direction perpendicular 
to the galactic plane) plays a dominant role 
during the early evolution of radio sources. 
The general form used for this gradient is 
p(z) = p (a /z)P ( [King 1972) ), where p 



W&^-ho are the reference density and the 



2.1. Evolution of radio doubles 

The dynamic evolution of extra-galactic 
double radio sources is characterized by 
observables such as the kinetic power 
provided to the lobes, the total extent 
of the source, the advancing velocity 
of the terminal hotspots and the outer 
lobes, and depends on the density gradi- 
ents of the ambient medium in the host 
galaxy along the path of the jets and 
the lobes. Modeling of the hotspots, the 
lobes, and the bow-shock has been done 
using self-similarity arguments to main- 
tain the overall shape of a symmetric dou- 
ble source (e.g., Begelman Cioffi 1989 



Begelman 19961 ICarvalho 19941 ICarvalho 1 998; 



IKaiser fc Alexander 1997] IKaiser fc Best 20071) . 

The radiative properties of the struc- 
tural components are determined by the 
balance between the particle energy den- 
sity, adiabatic losses resulting from expan- 
sion into lower pressures regions, and at 
later stages the evolving energy density 
of the magnetic fields in the hotspots and 
lobes, and the energy density of the Cosmic 
Microwave Background (CMB). The ra- 



scale length of the ISM and the IGM, re- 
spectively. Observational data of elliptical 
galaxies, galaxy groups and clusters sug- 
gest a double-^ model fIXue fc Wu 20001 
IFukazawa et al. 20041) with an innermost 
atmosphere of roughly constant density 
and a steepening density profile beyond 
ao- A transition at ao ~ 1 kpc from 
the galactic ISM (ft « 0) to the IGM 
(fi > 1.5) (IKaiser fc Best 2007j) is con- 
firmed by the observed discontinuity of 
the hotspot radius to arm length rela- 
tion (IKawakatu et al. 2009|) . However, ao 
can become significantly larger for certain 
galaxies. 

2.2. Stages of evolution of radio 
doubles 

The existing literature provides a fun- 
damental understanding about the dynam- 
ical evolution of FRII-like double radio 
sources. Although there is no full agree- 
ment about the actual parameterization 
of this evolution, four evolutionary radia- 
tive stages for a source with constant jet 
power P src may be identified that are re- 
lated to the dominant physical processes 
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(see IKaiser fc Best 20071) : 
Phase 1 - CSO - The radio luminosity 
evolution in this stage of CSOs and com- 
pact GPS sources is characterized by a 
flat density profile (up to a source size of 
about 1-3 kpc) with = 0. During this 
first stage adiabatic losses in the hotspots 
and lobes dominate the energy dissipation 
process as the embedded B-fields are not 
yet strong enough to cause dominant syn- 
chrotron radiation losses. The radio power 
increases with time as P ra( i oc t 2 / 5 and with 
the source size as P rad oc £) (8 ~ 7/3)/12 = D 2/3 
because of a steady increase in the con- 
version efficiency from jet kinetic power to 
radiative power. The spectrum should be 
steep with index a = 1.0. 
Phase 2 - MSO-1 - This stage of ex- 
tended GPS and more compact CSS sources 
occurs in the transition region (ao = 1 — 3 
kpc and /3 ~ 1) between the inner galaxy 
with a flat density profile and the outer 
galaxy with a steeper profile. During this 
phase there is a balance between adia- 
batic losses and synchrotron losses in the 
hotspots and lobes. The radio luminosity 
varies as P rad oc D°. The radio spectrum 
is again rather steep with a = 1.0. 
Phase 3 - MSO-2 - The third stage 
is mostly occupied with CSS sources for 
which the synchrotron losses resulting 
from entrained B-fields dominate the adi- 
abatic expansion losses resulting from the 
steep density gradient (a > 3 kpc and 
(3 > 1.5) of the external (inter-)galactic 
environment. The radio luminosity de- 
creases steadily with source size as P Ta( i oc 
£,(8-70/12 = £,-o.2 ; while the radio gpec _ 

trum would be relatively flat with a = 0.5. 
Phase 4 - LSO - The fourth evolution- 
ary stage is populated with fully developed 
FRII and FRI sources. The density pro- 



file of the IGM follows a power law with 
/3 > 1.5. The radio luminosity decreases 
sharply as P rad oc D^"^^ 5 -^ = D~ le 
and the radio spectra are steep with a = 
1.0. This phase may start around 100 kpc, 
where the inverse Compton losses resulting 
from the CMB energy density dominate 
over the synchrotron losses. 

The above described evolutionary stages 
have been depicted in the luminosity - 
size (Prad — D) diagram of Figure [I] for 
two representative power levels, that rep- 
resent a complete CSO-MSO-LSO lumi- 
nosity evolution. Along these evolution- 
ary tracks it is assumed that the jet power 
P src remains constant during the whole 
evolutionary lifetime until the LSO FRII 
stage is reached. The shape of this radio 
power P rad track is fully determined by the 
(/3 — D) relation. 

2.3. Jet stability and flaring sources 

The P ra d — D evolution of FR I-like flar- 
ing sources can partially coincide with 
that of the well-defined FRII-like dou- 
ble;sources. During the flaring stage the 
self-similar modeling for the morphology 
and the radiative properties would apply. 
Although both may coexist in the same re- 
gion of the (Prad — D) diagram, the (down- 
ward) evolutionary tracks of flaring sources 
are determined by the lower radio power 
and the cooling of the expanding radio 
structures. 

The stability of the jet flow and the in- 
teraction with the ambient material deter- 
mine the efficiency of the energy and mo- 
mentum transport to the lobes as well as 
the relative radio prominence of the jets 
and the lobes. After the oblique (reconven- 
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ing) shock at the base of the jet converts 
the expanding flow profile to a nearly par- 
allel (cylindrical) flow, a relativistic flow 
may persist over a long distance until even- 
tually a decelerating shock forms at the en- 
try of the lobe. High-power jets with lami- 
nar flow are most efficient for transporting 
the energy and momentum, which makes 
the jet mostly invisible in the radio maps 
of most prominent FR II-like doubles. 

A fluid boundary layer exists between 
the supersonic and relativistic core of 
the jet flow and the stationary ambient 
medium that forms the velocity gradi- 
ent from maximum to minimum veloc- 
ity. This boundary layer dissipates ki- 
netic energy because of viscous entrain- 
ment and fluid instabilities. With dis- 
tance traveled the boundary layer grows 
in thickness and at a certain characteristic 
scale length after the confining shock this 
initially laminar boundary layer becomes 
turbulent, which further increases the en- 
trainment and the momentum loss and 
weakens the jet. Eventually the bound- 
ary layer occupies the whole width of the 
jet. Under normal conditions, the su- 
personic jet with boundary layer forms 
a stagnation standing shock in the lobe 
flKaiser fc Alexander 19971 [Wang et al. 201 ip 



tance at which the boundary layer has 
grown to cover the whole width of the jet 
flKaiser fc Alexander 19971 [Wang fc Kaiser 2008 



which defines the location of the hotspot. 
However, a jet with a boundary layer 
reaching its centerline has lost much of its 
energy and momentum and fails to form a 
standing shock at a stagnation point. This 
jet ceases to supply sufficient energy and 
momentum to the lobe, which becomes 
diffuse and undefined as in FRI sources 



( ]Wang et al. 2009[ ). 

Stability analysis of the jet flow within 
the ambient medium predicts the dis- 



Wang et al. 2011| ). This distance may be 
expressed as: 

z e = Cl D^' G (1) 

where the proportionality constant varies 
with the jet power and the properties of the 
density profile (IKaiser Alexander 1997)) . 
This distance forms a specific boundary in 
the P ra d — D diagram, beyond which the 
jets become ineffective in supplying energy 
and momentum to the lobes. 

Noting that the shape of this boundary 
in the P ra d — D diagram (as a function of 
jet travel distance) remains model depen- 
dent, a representative boundary region has 
been entered in Figure [I] (denoted by two 
black dotted lines) using the same (3-D 
variation as used for the radio luminosity 
tracks. The vertical scaling of this region, 
which varies with the environmental pa- 
rameters and the jet power, is not well- 
determined and has been chosen to incor- 
porate some prominent FRI sources (3C31 
with z e = 3.5 kpc and P rad = 2.0e24 W/Hz 
and 3C449 with z e = 17 kpc and P ra d = 
1.12e24 W/Hz) and matches the FRII - 
FRI boundary at larger values of D. Evo- 
lutionary tracks reaching into this stabil- 
ity region predict jet stagnation/disruption 
and flaring FR I-like lobe structures. In ad- 
dition, this curve predicts that low-power 
CSO and MSO sources, or any source that 
experiences a systematic decrease (or ces- 
sation) of jet power, could directly evolve 
into a stagnated FRI-like source. 

2.4. Jet power continuity 

The continued evolution of radio sources 
described above depend primarily on the 
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continuity of the nuclear activity and the 
persistence of the jet. Theoretical stud- 
ies indicate that short-lived compact radio 
sources and the intermittent activity of the 
central engine may be caused by a radia- 
tion pressure instability within an accre- 
tion disk QCzerny et al. 2009 ). According 
to these authors, a radio source powered 
by a short-lived outburst of the central ac- 
tivity is not able to escape from the host 
galaxy unless the active phase lasts longer 
than 10 4 year, which is typically the time 
required for the lobe to successfully pass 
the ISM-IGM boundary (see discussion in 
sections 2.3 and 4.6). The discontinuity of 
jet activity has been associated with the 
formation of FRI-like sources from CSSs 
because of intermittent and episodic AGN 



activity (Marecki Spencer &: Kunert 2003 



Kunert-Bajraszewska et al. 2006[ ) and the 
death of CSO sources before reaching the 
MSO stage fISnellen et al. 1999( |Gugliucci efaT 



fuse, the hotspot separation will also be 
retarded and the source moves downward 
in the P ra d — D diagram. After re-ignition 
the retarded hotspot could again move to 
the entry of the lobe. A powerful recurrent 
source will expand its jet into the exca- 
vated channel of past activity and may ex- 
hibit a double-double morphology. Alter- 
natively, restart in the jet along a different 
axis may result in X-shaped source. A long 
time separation between subsequent events 
(compared with the cooling lifetime) will 
make the source reappear as a new startup 
FRII-like source. 

2.5. The size and velocity of the 
hotspots 

The size of the hotspot results from 
a balance between the energy carried by 
the jet and losses from radiation and adi- 
ibatic expansion. The parametric solu- 



Kunert-Bajraszewska et al. 2010) 



IBS* 



A reduction of the jet power results in 
the formation of a hotspot receding from 
the leading edge (and the entry point) 
of the lobe, the jet becomes more lossy 
because of increased prominence of the 
jet boundary, and the lobe develops an 
instability-driven (meandering) structure 
as in FRI galaxies. The emission spec- 
trum of the lobe will steepen and the jet 
may develop surface instabilities. Further 
reduction or termination of the jet power 
turns the source into a low-power relic as 
it eventually dies in its current phase. 

Recurrent jet activity will result in var- 
ious FRI-like morphologies depending on 
the jet power and the time separation be- 
tween subsequent events. During a low- 
power interval, the lobe becomes more dif- 



icris use a simplified assumption that 
the separation velocity of the hotspot re- 
mains constant during the early evolution 
fIFanti et al. 19951 IReadhead et al. 19961 
IGarvalho 19"98| IKaiser fc Alexander T997j) . 
However, constant-velocity models may 
not be representative for observed source 
behaviour. Analysis of the hotspot size 
versus D relation shows that the ob- 
served hotspot size for CSOs has a dif- 
ferent growth rate with arm length (as 
D lA for D < lkpc) than the larger MSO 
(GPS, CSS) and LSO-FRII sources (as 
D 0A for D > lkpc), which also confirms 
the location of the ISM— IGM transition 
(IKawakatu fc Kino 2006} IKawakatu et al. 2009|) 



Because the change in size of the hotspot 
and the hotspot separation velocity are 
parametrically related, these rates suggests 
that the hotspot of CSOs should decelerate 
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as Vhs D~ x (for (3 = 0) and accelerate 
as vhs °^ D 3 (for (3 = 1.5) during later 
phases. Under our simplified assumptions 
the vhs for the CSOs would decrease as 
D~ 2 / 3 during Phase 1 with constant P src 
(Section 4.1). The variation in the de- 



tected surface brightness levels for differ- 
ent observing arrays prevents a test of this 
relation using the CSO data in this paper. 

It should also be recognized that the 
jet remains at least supersonic in order 
to generate a hotspot. If vhs continues 
to decrease linearly until it reaches the 
ISM— IGM boundary, the initial hotspot 
advancing velocity should at least be 
0.3-0.5 c (see IKawakatu fc Kino 2006 



IKawakatu et al. 2009|) . Jet flows with ve- 
locities below this threshold become sub- 
sonic (of order 10~ 3 c) before reaching the 
ISM— IGM boundary and these will die 
prematurely (see section 4.4). 

2.6. Morphological Classification of 
Sources 

The physical processes that determine 
the evolution of radio sources create a large 
variety in the radio morphological struc- 
ture of radio sources. The following mor- 
phological types (Mtype) may be assigned 
to members of the family of extragalactic 
radio sources. 

Mtype 2 double sources have ener- 
getic laminar jets with prominent well- 
confined lobes but weak or invisible cen- 
tral core emission and weak or no interven- 
ing jet structure. They have well-confined 
lobes and prominent terminal hotspots at 
the leading edge of the lobe with a sharp 
edge-brightened morphology, that is in- 
dicative of a strong advancing shock. A low 
hotspot size— to— linear extent ratio indi- 



cates dominant radial expansion. Mtype 2 
sources have sufficient and long-lived jet 
power and occur among the CSO, MSO, 
and LSO populations. Mtype 2 sources 
evolve along the evolutionary tracks in the 
radio power-size (P ra d — D) diagram. FRII 
sources are the largest LSO sources within 
the Mtype 2 double population. 

Mtype 1 flaring sources are charac- 
terized by lossy prominent jets, diffuse, di- 
verging, and less confined lobes or flaring 
lobes, and hotspots located away from the 
lobes. Mtype 1 sources have jets where 
the boundary layer reaches the center- 
line (see section 2.3), which results in 



a low- momentum, flaring and/or mean- 
dering flow constituting the radio lobe. 
Mtype 1 sources exist among lower-power 
CSO, MSO and LSO sources that have 
reached the jet (in-)stability region in the 
P-D diagram (see Figure 1). Mtype 1 MSO 
and LSO sources are the FRI sources with 
hotspots at < 50% of the arm lengths. 

Mtype 3 dying sources experience 
cessation or reduction of the jet power, 
which result in the disappearance of the 
hotspots. The remaining radio structure 
continue its adiabatic expansion without 
the necessary energizing of the lobes. The 
sources becomes an expanding diffuse relic 
unless there is reactivation of the nu- 
clear activity and re-formation of the jet. 
Among CSOs the short-lived central source 
activity results in flow instabilities before 
the ISM— IGM transition point, and the 
eventual dissipation of the lobe structure. 
MSO and LSO sources will retain their 
shape without a hotspot and continue to 
expand while radiating their energy away. 

Mtype 4 restart sources experience a 
restart of nuclear activity or intermittent 
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activity and the startup of jet power into 
the earlier excavated channels. Depend- 
ing on the new jet power, the hotspot ad- 
vances through the channel until it reaches 
the lobe entry and sets up a new shock- 
confined lobe within a larger existing lobe. 
Depending on the time interval between 
the periods of activity, the new lobe may be 
quite distinct from the relic lobe structure, 
as in double-double radio sources. Regu- 
lar intermittent activity may results in a 
series of shock fronts passing through the 
lobe. X-shaped radio sources could also 
signify a restart of the central source with 
an different jet axis. 

Mtype 5 obstructed sources have 
lower-power jets in the CSO Phase I and 
MSO Phase II that interact with the sur- 
rounding medium along its path, which re- 
sults in hotspots and misaligned outflows 
(ILonsdale &: Barthel 19861) . in analogy to 
the dentist's drill model (IScheuer 1982ft or 
the frustration model, and significant side- 
ways motions dissipating kinetic energy. 
A higher hotspot size— to— linear extent 
ratio implies significant sideways expan- 
sion/motions. The morphology of these 
sources will deviate from the straight two- 
lobe picture to be expected from Mtype 2. 

Mtype 6 core-jet sources have a 
small angle between the jet-direction and 
the line of sight. These sources appear as 
one-sided jet with enhanced radio bright- 
ness because of doppler boosting. 

Mtype 7 trail sources are trail radio 
galaxies where the Mtype 1 lobe structures 
outside the ISM— IGM boundary are bent 
because of ram pressure resulting from the 
motion of the host galaxy throughout the 
intra-cluster medium. 



2.7. Evolutionary scenarios 

The long-term evolutionary tracks through 
the P ra d — D diagram depicted in Figure [T] 
reflect the long-term Mtype 2 CSO-MSO- 
LSO double source evolution. All dou- 
ble radio sources will start out on such a 
track and they will continue to follow it as 
long as their jet power remains at the ap- 
propriate level. These evolutionary tracks 
also indicate that the Mtype 2 FR II double 
sources will naturally evolve into Mtype 1 
flaring FRI sources, because they eventu- 
ally reach the region of jet instability (see 



Section 2.3). This evolutionary path from 



FR II to FR I has been depicted graphically 
(blue arrows) in Figure [2j 

Cessation or changes of the jet power 
could happen at any moment during the 
CSO-MSO-LSO sequence. Mtype 3 dying 
sources undergoing a lowering or cessa- 
tion of the jet power move downward (red 
dashed arrows) in the P ra d — D diagram 
and may be found among the Mtype 2 dou- 
ble sources. Mtype 4 restart sources ex- 
periencing an increase in power will move 
upward in the diagram (black dashed ar- 
rows) to a higher power Mtype 2 source. 
Re-energizing a Mtype 3 source may again 
turn the source into a Mtype 2 double. 

Mtype 5 obstructed sources with some 
signs of disruption may be found through- 
out the whole CSO-MSO region in the 
P rad - D diagram. 

All sources having reached the jet stabil- 
ity region are already destined to become 
Mtype 1. 

The scenarios for radio source evolution 
presented here are distinctly different from 



those presented by (Kunert-Bajraszewska et al. 2010) 



These authors focus on the evolution from 
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CSS to FRII and FRI and do not address 
the earliest evolutionary CSO stages con- 
sidered in this paper. While the long-term 
evolutionary scenario for Mtype2 sources 
is similar for the scenario in Kunert- 
Bajraszewska et al. (2010), the details and 
the origin of Mtype 1 sources is different. 
Because Mtype 3 sources occur naturally in 
the whole diagram, the formation of low- 
power MSO-CSS sources does not (only) 
have to result from multiple interactions 
of a high-power CSS with obstacles in the 
host galaxy. Intrinsically low-power CSOs 
naturally evolve into low-power CSSs and 
also a reduction of the jet power in a CSO 
will result in a low-power CSS. 

3. The CSO data sample 

In order to compose a more complete 
sample of the CSOs for the purpose of a 
kinematics study, observations have been 
made of a sample of 10 CSO candidates 
without reliable proper motion measure- 
ments. The new observations were car- 
ried out at five frequency bands ranging 
from 1.6 to 15 GHz with the Very Long 
Baseline Array in 2005, and with a com- 
bined array of the Chinese VLBI Network 
(CVN) fILi et al. 20101) and the European 
VLBI Network (EVN) in 2009. The newly 
obtained data combined with the archival 
VLBI data made it possible to determine 
the separation velocities of the hotspots 
or to set upper limits for all ten sources 
(lAn et al. 2012p . Six of these sources are 
identified as CSOs and one as a CSO can- 
didate, which increases the sample of CSOs 
and candidates with known proper motions 
from 30 to 37. 

The classification of CSOs is based on 
the morphology and spectral index cri- 



teria. A CSO is defined clS 8b compact 
source (overall size <1 kpc) with two ex- 
tended emission components on either side 
of a central flat-spectrum core. While the 
central core may be too weak for VLBI 
imaging, the presence of two extended 
components with steep spectra and edge- 
brightened morphology would be enough 
for CSO classification. The mirror symme- 
try of the proper motions of the hotspots 
and internal jet knots with respect to the 
geometric centre provides supplementary 
CSO identification. Core-jet sources have 
not been included in the sample, because 
their radio luminosity is affected by an (un- 
known) Doppler boosting factor. A known 
redshift is required to convert from pro- 
jected angular size to physical linear size. 

A sample of 46 CSOs with known red- 
shift and size information has been pre- 
sented in Table 1. A total of 24 sources 
in this sample also have published CSO 
proper motion data for the purpose of a 
statistical study of the kinematic proper- 
ties of CSOs. Table 1 contains the follow- 
ing entries: 

Column (1) gives the object name. 

Columns (2) and (3) present the red- 
shift z and luminosity distance Dl, which 
are taken from the NASA Extragalactic 
Database (NED) and from the literatures. 

Columns (4) and (5) present the overall 
source size in both angular size (9as) an d 
projected linear size (LS). The calculation 
of Dl and the conversion from 9 as to LS 
make use of a cosmological model with Ho 
= 73 km s" 1 Mpc" 1 , tt M = 0.27, and H\ 
= 0.73. 

Column (6) gives the angular separation 
velocity of the hotspot n = 9 as I 'Ai, which 
is calculated as the rate of the separation 
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of two hotspots over the time interval. For 
most sources, fi is determined as a rela- 
tive separation velocity with one hotspot 
as the reference. For CSOs with a visi- 
ble core, the proper motions of hotspots 
are calculated as the separation velocity 
of each hotspot with respect to the cen- 
tral core. When multiple measurements 
are available, the value with the least un- 
certainty has been used. 

Column (7) gives the apparent trans- 
verse velocity in the source rest frame 
v in units of c obtained from /i: v = 
0.0158^/^(1 + z), where Da is the an- 
gular size distance of the source in Mpc 
and \x is the angular separation velocity in 
mas yr _1 . 

Column (8) presents the kinematical age 
in the source rest-frame, 9as/^(1 + z) or 
LS/v. 

Column (9) and (10) list the morphol- 
ogy classification based on the VLBI im- 
ages and the reference for each object (see 
details of the morphology classification in 
Section 2.5). 

Columns (11) to (13) present the flux 
densities S at 1.4, 4.8 and 8.4 GHz, re- 
spectively. 

Column (14) presents the calculated 
turnover frequency. 

Column (15) presents the spectral in- 
dex a|;| between 4.8 and 8.4 GHz (defined 
as S u oc z/ _Q ). A sample of seven CSOs 
were monitored using the VLA at 8.5 GHz 
and they showed little variability (mean 
rms of 0.7%) over a period of eight months 
dFassnacht fc Taylor 200lj ). Other CSOs, 
such as OQ 208, display much larger flux 



frequency flux density measurements of 
CSOs have been used preferentially from 
two available datasets: VLA A-array data 
at 0.3, 1.4, 4.8 and 8.4 GHz at three epochs 



during 1984 and 1991 dStanghellini et al. 1998Q , 
and VLA data at six frequency bands 
ranging from 1.4 to 43 GHz in 2003-2004 
flOrienti et al. 20071) . For CSOs without 
simultaneous data, flux densities have been 
used from the Green Bank Telescope at 
1.4 and 4.85 GHz flWhite fc Becker 19921 
Gregory fc Condon 1991]) and the VLA-A 



at 8.4 GHz flPatnaik et al. 19921 |Healey et al. 2007[ ) 

Column (16) lists the absorption-corrected 
radio power at 1.4 GHz in units of W 
Hz^ 1 determined as Pxaghz = 4rD\(1 + 
z)^ 1 ' Siaghz- CSOs often exhibit in- 
verted spectra with a turnover at a few 
GHz dStanghellini et al. 1998[IOrienti et al. 2007p 



caused by synchrotron self-absorption and/or 
free-free absorption (FFA). The observed 
flux densities in Columns (11)- (13) show 
an Siaghz much lower than the 1.4-GHz 
flux density extrapolated from S^ghz 
and the spectral index afg, which is in 
agreement with a spectrum turnover at 
frequencies higher than 1.4 GHz (Col- 
umn 14). Since can be a reason- 
able representation of the optically-thin 
section of the spectrum, the 4.8-GHz 
flux density (S^sghz) and spectral index 



a 



8.-1 



have been used to extrapolate the 
absorption-corrected 1.4-GHz flux den- 
sity used for computing P\aghz- Those 
sources with a turnover much higher than 
5 GHz, the 8.4-GHz flux density and the 
spectral index between 8.4 and 15 GHz 



have been used ( Stanghellini et al. 1998 



variations(de Bruyn 1990 Stanghellini et al. 1997 



IQrienti et al. 2007j) . 



IWu et al. 20121) . Quasi-simultaneous multiple- 
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4. Dynamic Evolution of CSOs 

The smallest and youngest members of 
the radio source population are the CSO 
sources with linear sizes less than 1 kpc. 
All radio source must go through this stage 
where they are still submerged in the dense 
nuclear ISM of the host galaxy, which 
provides the strongest resistance against 
growth. The CSO sources occupy the left 
side of the radio power-size P ra d ~ D dia- 
gram (Fig. 1) and show a range of five or- 
ders of magnitude in radio power. Because 
the adiabatic losses dominate over the syn- 
chrotron losses for these compact sources, 
the P src remains larger than P ra d until the 
ISM— IGM boundary radius a^. The mor- 
phological type to be expected during this 
early evolutionary stage ranges between 
Mtype 2 doubles, Mtype 1 flaring sources 
at the lower end of the P ra( r range, short- 
lived Mtype 3 dying sources, and Mtype 5 
obstructed sources. The dynamical proper- 
ties of the CSO sample will thus provide 
essential information about the evolution- 
ary stage of the sources. 

The observable parameters that may 
be used to study the dynamic properties 
of the CSO sources are (1) the hotspot 
separation velocity (vhs equals twice the 
hotspot advance velocity) assuming the 
jets are ejected with the same velocity on 
opposite sides, (2) the overall linear size of 
the source (D), and (3) the radio power 
(Prad), which will equal the total jet power 
P src at the ISM— IGM transition distance. 
The dynamic behavior of the CSO sources 
are subject to variation of the jet kinematic 
and radiative powers, the variation in the 
ISM and IGM density gradients, and the 
initial density as well as scale height. Pro- 
jection of the source structure on the plane 



of the sky will lower the observed values of 
the projected source size and the hotspot 
separation velocity. The radio power and 
the kinematic age of the source are not af- 
fected by projection. 

4.1. Modeling the CSO stage 

The evolution of the lobes and the 
hotspots of sources during the CSO stage 
have been modeled assuming self-similarity 
and expansion into a fixed opening angle. 
Assuming that P src is constant during the 
CSO phase, this results in an overall mo- 
mentum balance equation at the bow shock 
flKaiser fc Alexander 1997) : 

P src = c 2 D 2 -?v% s (2) 
where P src is the source mechanical power 
and vhs is the hotspot separation velocity, 
assuming that the this velocity varies pro- 
portionally with the bow shock and lobe 
separation velocities. In addition, the lin- 
ear extent D varies with evolutionary time 
t as flKaiser fc Best 20071) : 

D = c 3 (P src /(pa^)r-H 3 /^ (3) 

During the CSO stage, the radio power is 

less than P src and varies as (IKaiser fc Best 2 007): 

The proportionality constants in these 
equations depend on the properties of the 
ISM and the power of the source as pre- 
sented in the above references. 

The above equations with /3 = may be 
used to describe the dependences of the ob- 
servables presented in this paper, these be- 
ing the distance of the lobe D, the hotspot 
(or lobe) separation velocity vhs, the ob- 
served radio power P ra d, and the evolution- 
ary time t measured as 'kinematic age" . 
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Fig. 3. — Variation of the CSO radio power 
with projected linear source size. Symbols 
indicate high (red circles) and low (black 
squares) radio powers. Sources with con- 
stant jet power P src have a predicted evo- 
lution P ra d oc D 2 / 3 indicated by the dashed 
arrows. The general trend in the distribu- 
tion is represented by a P ra d oc D A relation 
(dotted line). 

4.2. Radio Power versus Linear Size 

More powerful sources should be larger 
for a certain evolutionary age. However, 
there is no uniform relation between radio 
power and so urce size. The modeling re- 
sults (Section 4.1) suggests a P ra d oc D 2 ^ 3 
relation depicted as dashed arrows in Fig- 
ure [U This relation reflects the increase 
of the conversion efficiency from the ki- 
netic power of the jet to the radio power as 
sources get larger and synchrotron losses 
increase. The spread in the data follows 
from variation of the source power P src and 
the ISM density structure. Projection ef- 
fects would shrink the apparent source size 
and move a source to the left of the dia- 
gram. 

More powerful sources can indeed be 
larger but they may also be older. The 
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Fig. 4. — Variation of the CSO radio power 
versus the hotspot separation velocity. As 
a group, the sources with higher vhs gen- 
erally have a higher radio power P ra d and 
also a higher source power P src - Individual 
sources are predicted to evolve from high- 
ly and \ow-P rad to low-t> HS and high-P md 
with time as P ra d oc vjl s . Symbols in- 
dicate high-power (circles) and low-power 
(squares) CSOs. Arrows on the symbols 
indicate upper limits. 

relation between the source size and kine- 
matic age in Figure [^-b confirms that older 
sources are in general larger in size. The 
lower power sources (squares) reveal this 
trend in Figure [3] High-power sources (cir- 
cles) show a different distribution, possi- 
bly resulting from source selection effects, 
but they may still behave according to the 
same relation. 

The radio power is the measured quan- 
tity, while the physics of the source is de- 
termined by the source power P src . Be- 
cause P ra d would equal P src at distance 
D = cio = 1 kpc, an extrapolation using 
P ra d oc D 2 ! 3 to this point will provide an 
estimate for the source power P src . 
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4.3. Radio Power versus Hotspot 
Velocity 

More powerful sources should have 
larger hotspot separation velocities for a 
given environment and a given age. How- 
ever, the diverse environments and ages of 
the CSO sample cause a spread in veloc- 
ity in Figure |4j A simple relation follow- 
ing from momentum conservation and the 
variation of P ra d with D (Sec. 4.1) sug- 
gests that a source with fixed P src varies 
as P rac [ oc v]j S , i.e. a source with a given 
P src will evolve to larger P Ta d and smaller 
vhs until it reaches the ISM— IGM tran- 
sition point. As a group, the CSOs with 
higher vhs a ls° have higher P ra d, because 
of their higher source powers P src - 

4.4. Hotspot Velocity versus Linear 
Size 

The dynamics of the evolution of the 
radio source within the varying ambient 
environment should be well depicted by 
the variation of hotspot separation veloc- 
ity with projected linear size. The vhs — D 
variation should be indicative of the valid- 
ity of the parametric similarity solutions. 

Modeling predicts that the hotspot sep- 
aration velocity should vary as vhs oc 
£)- 2 /3 f or a source with constant source 
power P src (Fig. [5]). The separation veloc- 
ity of an individual source expectedly de- 
creases as the source becomes larger. The 
D^ 1 dependence derived from the varia- 
tion of the hotspot size is also consistent 
with these data (IKawakatu fc Kino 2 006; 
IKawakatu et al. 20091) . As a group, the 
more powerful sources with faster moving 
hotspots will also be larger. This property 
may be represented by a vhs oc D 2 ^ de- 
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Fig. 5. — Variation of the hotspot sepa- 
ration velocity versus the projected linear 
size. The dashed arrows vhs °c D~ 2 ^ in- 
dicate the predicted evolution for a source 
with constant power. The dotted line rep- 
resents the lower boundary for the jet flow 
to continue to accelerate at the ISM— IGM 
transition point. CSO close to and below 
this line will not survive as double sources 
beyond the ISM. Projection effects for the 
CSOs shift data points diagonally towards 
the origin. 



pendence. 



for a ra- 



As discussed in Section [2T3 
dio source evolving along the primary 
evolution track from CSO to CSS and 
to FRII stage, the hotspot separation 
velocity must at least be supersonic at 
the transition ISM— IGM transition point 
(IKawakatu fc Kino 2006| IKawakatu et al. 2009ft 



This means that a (survivor) CSO should 
lie above the (dotted line) vhs oc D~ 2 ^ 3 
threshold in Figure [5j CSO sources cur- 
rently below this threshold are Mtype 1 
flaring or Mtype 3 dying sources. 
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Fig. 6. — Variation of source radio power 
with kinematic age. Modeling predicts 
an evolutionary relation of P ra d oc 
(dashed arrows). 



4.5. Variation with Kinematic Age 

The hotspot separation velocity, the lin- 
ear extent, and the radio power of the 
source all vary strongly with kinematic age 
(or evolutionary time). Modeling results 
assuming a constant jet power P src during 
Phase 1 predict an increase of the radio 
power with evolutionary time P ra d oc T%(% 



(Sect. 4.1). This relation may explain 
the data points using the grouping in the 
power levels (Fig. [6]). 

The hotspot velocity shows a large 
spread and are consistent with a decrease 
with Kinetic Age in Figure [7k. Modeling 

2/5 ^""^ 

predicts a vhs °c T kin ' evolution for a 
source with fixed power P src and decelera- 
tion for the hotspot with time for CSOs. 

The linear size of a CSO naturally in- 
creases with the Kinematic Age as is man- 
ifested in Figure [7)d. Modeling predicts 
that D oc T^ 5 for an individual source 
with constant P src (dashed arrows) (see 
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Sect. 4.1). The general trend that older 
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Fig. 7. — Variation of the linear size and 
the separation velocity with kinematic age. 
(a) Both source groups (high and low ra- 
dio power) are consistent with hotspot sep- 
aration velocity decreasing with age. The 
predicted evolution of the sources follows a 
Vhs oc T ki ^ 5 relation (dashed arrows), (b) 
The CSO linear size shows an expected in- 
crease with age. The predicted variation 
for individual sources with constant power 
is D oc as indicated with dashed ar- 
rows. As a group also the older sources 
have larger sizes as D oc (dotted line). 



sources have larger sizes is manifested by 
the whole group and may be expressed 
with a D oc Tjjj^j as indicated with the dot- 
ted line. 
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4.6. Spectral Index versus Time 
and Size 

The spectral index will vary strongly 
when a CSO evolves from (smaller) adi- 
abatic expansion dominated sources to 
(larger) synchrotron-loss dominated sources. 
Synchrotron self-absorption will dominate 
for these small-size CSO sources. Mod- 
eling of the radio properties suggests that 
the spectral index for CSOs should be close 



to unity (Sect. |2.2). 



In the diagrams of Figure [8| the spectral 
index a|;|, defined as S u oc v~ a , is plot- 
ted against kinematic age T kin and source 
size D. The spectrum becomes steeper for 
larger kinematic ages and linear sizes as the 
synchrotron self-absorption opacity gradu- 
ally decreases. Both diagrams display this 
expected general increase from a = 0.3 to 
unity during the CSO stage. Sources that 
do not follow the trend include a single 
(relic) source above the distribution and 
four High Frequency Peakers with a very 
high turnover frequency (y to >5 GHz) be- 
low the distribution. 

The spectral index for the CSO sample 
presented here is taken in the optically- 
thin part of the spectrum above the 
turnover frequency of the sources, which 
decreases systematically with sources size 
(IQ'Dea 1998|) . The observed range of of 
spectral indices above the peak for samples 
of CSS and GPS sources is a = 0.5 — 0.9 
( IQ'Dea et al. 19 90). which partially over- 
laps with the CSO range. The same 
range was derived for the spectral in- 
dex above spectra peak for GPS and FR 
galaxies, suggesting that particle accel- 
eration and energy loss mechanism pre- 
serve the same average spectral index 
over most of the lifetime of the source 



fide Vries et al. 2009|) . Variability of the 
spectral index has also been associated 
with the HFP sources and to identify truly 
young sources (IQrienti et al. 2007j) . 

The observed evolution of the spectral 
index of the CSO sample confirms the (ex- 
pected) systematic change as the source 
gets older and larger, which is explained 
by the shift of self-absorption to lower fre- 
quencies. The youngest CSO and HFP 
sources have the flattest spectrum. The 
observed CSO trend should continue into 
the MSO stage, where it should level off to 
continue with a (a = 1) steep spectrum. 

5. Statistics of the CSO sample 

This unbiased sample of CSO sources 
used in this paper is found not to be uni- 
formly distributed in power, age and size. 
In general, the CSO phase represents a 
short episode in the life of a radio source 
and a uniform distribution of CSO sources 
of all ages would produce a flat distribution 
with age. However, the histogram of the 
Kinematic Age shows a steadily decreas- 
ing distribution (Fig. [9j), which suggests 
that some CSOs will not become old. The 
halfwidth of this distribution is on the or- 
der of 1100 years. The relative excess of 
young sources confirms the earlier results 
of ( JGugliucci et al. 2005] lAn et al. 20121) . 

Similarly, the histogram of D shows a 
decreasing number of sources having larger 



sizes (Fig. 10). Projection effects can not 
be the cause of this distribution because 
for a randomly oriented sample, the prob- 
ability for a source to be close to the plane 
of the sky is much larger than for it hav- 
ing a large inclination angle (i.e, pointing 
at the observer). The probability distribu- 
tion for a distribution with random inclina- 
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tion angles and a given size D increases lin- 
early until it reaches D (a triangular distri- 
bution). Similarly, a uniform distribution 
with random source sizes would also give a 
flat distribution if all grow larger into CSS 
sources. The observed D distribution indi- 
cates a halfwidth of about 100 pc for the D 
histogram. The flat distributions at higher 
D and Age suggest that sources older than 
typically 1100 years and larger than 100 pc 
have a better survival rate and will evolve 
beyond the CSO stage. 

The excess of sources at small sizes and 
Kinematic Ages indicate that a significant 
fraction of CSO sources does not continue 
to evolve beyond the CSO stage towards 
the MSO stage. This depletion of CSOs 
may result from reduction/cessation of the 
jet power and the ensuing instabilities, or 
from total termination of the nuclear ac- 
tivity after a short activity period. In ad- 
dition, any recurrent or periodic ejection 
of emission components would result in 
the re-excitation of an existing lobe cavity 
created during earlier outbursts but even 
such sources may not grow into Mpc-scale 
doubles. This implication is supported by 
CSO observations showing sub-pc-scale ex- 
tended relics, for example, 1511+0518 and 
OQ 208 (? IWu et al. 2012]) . 

5.1. Transient CSO sources 

Inspection of the source morphology us- 
ing the above considerations of a larger 
sample of 46 CSO sources with high- 
dynamic-range VLBI images from the lit- 
erature suggests that only 13 (29%) sources 
are Mtype2 double sources with well- 
confined hotspots at the leading edge of 
the lobes. Similarly within the current 
sample of 24 sources with kinematic data, 



a total of 7 (also 29%) sources have such 
a well-defined Mtype2 signature. The re- 
maining sources show decaying structures 
with prominent jets, hotspots before the 
lobes, and diffuse or disrupted lobes. Sev- 
eral sources may even qualify as Mtype3 
radio relics. 

In principle, the evolutionary character- 
istics provided in this paper would also al- 
low the identification of the evolutionary 
state of the CSOs. For instance, the largest 
sources with the lowest expansion veloci- 
ties could be precursors of Mtype ,1 flaring 
or Mtype 3 dying sources (in Fig. [5]) and a 
high separation velocity at a high age may 
indicate survival of the CSO stage. How- 
ever, the sources in our sample classified as 
Mtype 2 doubles do not present any clear 
tendencies or systematics in the evolution- 
ary diagrams of this paper. Therefore, 
the current evolutionary data do not (yet) 
provide discriminatory guidance about the 
chance of survival of an individual CSO 
and its morphological classification. 

5.2. Transient behaviour 

The long-term AGN duty cycle of 10 7-8 
yr is generally thought to be induced by 
major mergers and long-term accretion 
events. However, the current data suggests 
that a significant fraction of the CSO pop- 
ulation consists of short-lived radio sources 
with a typical lifetime less than 1000 years. 
Similarly, histograms of GPS and CSS 
source sizes show flat distributions fol- 
lowed by a steep distribution for FRI or 
II sources (IQ'Dea 19980 . Therefore, some 
GPS and CSS sources will not evolve into 
classical doubles and instead become tran- 
sient or frustrated f lO'Dea &: Baum 1997j) . 
A strong similarity may thus exist with 
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X-ray/optical/radio transient events in 
AGN resulting from tidal disruption of 
stars close to the supermassive black hole 
(SMBH). 

Short-term AGN duty cycles of < 10 4 yr 
may result from transient nuclear accretion 
events or instabilities in the accretion disk. 
CSO radio power of ltf 7 W/Hz (at 5 GHz) 
can be sustained by SMBH accretion when 
the accreted mass is: 

AM = 1.5 x 10- 3 (ar7) -1 AT M , 
where AT is the time interval in years, 
a Ri 0.1 represents the accretion efficiency 
with half the available mass accreting, and 
7] m 0.01 represents the overall conversion 
efficiency from accreted energy to radiative 
energy. A medium power CSO would thus 
require about 150 M for a 100 year tran- 
sient episode. Such numbers make source 
transience more plausible and may account 
for the disappearance of young CSOs. 

6. Conclusions 

In order to classify the energetics and 
morphology of CSOs and other radio 
sources, seven morphological types may 
be identified that also incorporate the 
large-scale FR I and FR II phenomenology. 
Lower power Mtype 1 flaring sources are 
FRl-like sources that show symptoms of 
jet decay, loss of power, diffuse lobes, or 
large-scale jet instability. They exist below 
the transition region in the P ra d ~ D dia- 
gram, unless they are re-energized. Higher 
power Mtype 2 double sources are FRII- 
like sources with efficient (often invisible) 
jets and well-confined hotspots and lobes. 
Such sources will follow the complete evo- 
lutionary path in the P ra( i — D diagram 
until higher D. Once an Mtype 2 source 



reaches a stage where the jets become un- 
stable, either by growth of instabilities or 
because of insufficient power, this Mtype 2 
source turns into an Mtype 1 source. 

Cessation or intermittency of jet power 
can happen during any evolutionary stage, 
and a Mtype 3 dying source will move 
downward in the P ra d ~ D diagram (see 
Fig. [2|. Therefore, throughout the whole 
Prad—D diagram there may be Mtype 3 dy- 
ing sources as relics of burned-out sources. 

Mtype 4 restart sources re-energize Mtype 3 
dying sources or even low-power Mtype 1 
flaring sources. The morphology of the 
Mtype 4 restart sources will depend strongly 
on the state of the sources before the 
startup. The morphology will differ strongly 
if the restart occurs in an Mtype 2 double 
source or in an Mtype 3 relic source. X- 
shaped radio sources would also fit in this 
morphological type. Mtype 5 obstructed 
sources occur among the CSO and MSO 
population and show evidence of blockage 
or bending of the jet flow, an example is the 
CSS 3C 48 flAn et al. 2010j) . Mtype 6 core- 
jet sources are subject to strong projection 
effects, while host galaxies of Mtype 1 trail 
sources have a significant motion through 
the intra-cluster medium. 

CSOs show a compact-double radio 
structure that is analogous to classical ex- 
tended doubles, although their sizes can 
be as much as five magnitudes smaller. 
The morphological similarity among these 
radio sources invokes an evolutionary sce- 
nario from CSO via MSO to LSO. Their 
evolution is characterized by the relations 
between radio luminosity, the separation 
distance of the two hotspots, the advanc- 
ing velocity of the hotspots, and the age 
of the source. The observable properties 
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of small-scale CSOs illustrate the factors 
that govern their dynamical evolution and 
eventually their fate, and the information 
of the physical properties of the ISM of 
the host galaxies on sub-kpc scales. For 
the purpose of a dynamic evolution study, 
an unbiased but incomplete sample of 24 
CSO sources has been selected solely on 
the basis of the availability of redshifts 
and hotspot separation velocities. 

Self-similar evolution models have been 
used to depict the dynamic evolution of 
extragalactic Mtype 2 radio sources in four 
stages i.e., CSO-GPS-CSS-FRII. The 
CSOs are in the earliest growing stage of 
radio sources where the radio power varies 
strongly with source size and adiabatic ex- 
pansion dominates over synchrotron losses 
inside the ISM-IGM transition point where 
the radial dependence (/?) of the ISM den- 
sity changes. This distance falls in the 
range of 1—3 kpc but it can also be sig- 
nificantly larger. The initial CSO sample 
indeed shows systematic behavior that is 
consistent with self-similar model predic- 
tions (e.g.. Kaiser fc Best 20071 

In all diagnostic diagrams using the 
CSO observables, the evolutionary path 
of individual sources can be parameterized 
by the predictions of self-similar modeling, 
assuming that the source maintains a con- 
stant jet power. According to these models 
an individual CSO increases its luminosity 
with evolutionary time and with size until 
it reaches the ISM— IGM transition region. 
Individual sources have a decreasing sepa- 
ration velocity with increasing power, with 
evolutionary time, and with size until the 
ISM— IGM transition. The sources system- 
atically increase in size with evolutionary 
time. 



The spectral index of the CSO sources 
changes from zero (flat) to unity (steep) 
during the course of their evolution until 
the ISM— IGM transition. This aspect is 
related to the self-absorption during the 
earliest stages of evolutions and is not cov- 
ered by self- similarity models. 

The CSO sample whole is not well- 
confined, because sources with different jet 
powers and different environmental param- 
eters are superposed in the diagrams. Nev- 
ertheless, the sample whole shows 
global variation that differs from the evo- 
lutionary path of individual sources. More 
powerful sources are larger and have higher 
expansion velocities, but cover the same 
range of kinetic age. Sources with higher 
expansion velocities are larger and larger 
sources are also older. 

In these models the hotspot separation 
velocity of a CSO with constant jet power 
varies as D~ 2/3 within the (3 = region. 
Beyond the ISM-IGM transition the veloc- 
ity in the jet may increase again if the jet 
flow remains laminar and supersonic. A 
sources that do not fulfill these conditions 
will not have a sustained growth beyond 
the ISM-IGM transition. 

The histogram of the CSO age esti- 
mates shows a lack of old sources, indi- 
cating that a significant fraction of CSOs 
will not evolve into MSOs and eventually 
LSOs. Sources with current age of > 1 100 
years and size of > 100 pc would have a 
better chance of survival. Classification 
of the current larger CSO sample sample 
shows that only 29% of the sample have 
clear Mtype 2 double qualities and these 
are likely to grow further in size. 

The dynamic properties and evolution- 
ary behavior of CSOs represent both the 
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youth and frustration scenarios. CSOs 
are naturally young radio sources at their 
early evolutionary stages. The CSOs with 
short-lived or intermittent nuclear activ- 
ity and/or having lower-power can be re- 
garded as frustrated or obstructed and will 
not evolve into full-size Mtype2 double 
sources. 

Further study of the dynamics of CSO 
requires a more complete sample of sources, 
because of the observational bias against 
the numerous high frequency peakers at 
the earliest CSO stages, particularly at 
higher power (compare Fig. j3J). The most 
compact radio sources are only identified 
as CSOs when the radio structure is re- 
solved at the VLBI resolution. As a re- 
sult, only sources with a hint of mas-scale 
structure from earlier VLBI surveys have 
warranted follow-up observations. A sys- 
tematic unbiased VLBI survey of a com- 
plete GPS sample would thus be necessary 
to identify the most compact CSO popula- 
tion. In addition, current CSO kinematic 
samples focus more on the high-power pop- 
ulation and more sensitive VLBI observa- 
tions of low-power CSOs are required for a 
complete view of the CSO evolution. 
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Table 1 

Radio Properties of a sample of 46 CSOs and Candidates 



Source 


z D L 


e AS 


LS 




V 


Age 


ID 


Ref. 


SlA 


•54.8 


Ss.4 


Vto 


a 8.4 
Q 4.8 


logioPi.4 


Ref. 




(Gpc) 


(mas) 


(pc) 


(/uas/yr) 


(c) 


(yr) 


(°) 


(Map) 


(Jy) 


(Jy) 


(Jy) 


(GHz) 


(") 


(W/Hz) 




(1) 


(2) (3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ii) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 



0003+2129 


0.455 


2.519 


3.8 


21.9 










O06 


0.100 





.260 


0. 


227 


5.4 


0.23(0.9) 


26.92 


O07 


0005+0524 


1.887 


14.641 


1.7 


14.5 










O06 


0.166 





.210 


0. 


.165 


3.7 


0.43(0.6) 


27.91 


O07 


0029+3456 


0.517 


2.854 


32.0 


200.0 










K04 


1.750 


1 


.312 


0. 


981 


0.4-1.4 


0.52 


27.02 


W92,G91,H07 


0038+2302 


0.096 


0.419 


18.5 


31.3 


41.0 


0.25 


410 


CI 1 


P09 


0.532 





.247 






0.4-1.4 


0.62 


25.03 


W92,G91 


0048+3157 


0.015 


0.0582 


4.4 


1.2 










VCS 


0.270 





.802 


0. 


286 


1.4-5.0 


1.88 


24.53 


W92,G06,H07 


0111+3906 


0.668 


3.917 


5.8 


39.0 


9.3 


0.34 


380 


CI 1 


OC98 


0.509 


1 


.301 


0. 


.918 


4.7 


0.62(1.1) 


28.10 


O07 


0119+3210 


0.060 


0.255 


89.4 


98.0 


<170.0 


<0.64 


>500 


CI I 


G03 


2.826 


1 


.571 


1. 


087 


0.4 


0.67 


25.44 


W92,G91,H07 


0428+3259 


0.479 


2.680 


2.7 


16.0 










O06 


0.148 





.524 


0. 


531 


6.9 


-0.02(0.6) 


27.06 


O07 


0431+2037 


0.219 


1.045 


42.7 


145.6 










F00 


3.611 


2 


.300 


1. 


.710 


0.63 


0.54 


26.73 


W92,G91,H07 


0503+0203 


0.584 


3.322 


10.36 


66.5 










S01 


2.118 


2 


.197 


1. 


411 


2.5 


0.81 


27.59 


W92,G91,H07 


0650+6001 


0.455 


2.519 


3.2 


18.5 










O06 


0.507 


1 


.150 


0. 


975 


5.3 


0.29(0.5) 


27.18 


O07 


0713+4349 


0.518 


2.870 


24.3 


146.0 


17.0 


0.51 


930 


CI 1 


PC03 


1.940 


1 


.670 


1. 


280 


1.9 


0.48 


27.38 


S98 


0943+5113 


0.42 


2.289 


3.9 


21.5 








CI II 


OD12 


0.077 





.147 


0. 


064 


3.7 


1.48(1.8) 


26.10 


OD12,ODS10 


0951+3451 


0.29 


1.482 


4.8 


20.7 








CI II 


OD12 


0.019 





.062 


0. 


055 


6.0 


0.21(0.6) 


25.66 


OD12,ODS10 


1035+5628 


0.450 


2.483 


32.2 


181.8 


37.6 


0.96 


620 


CI I 


T96 


1.780 


1 


.270 


0. 


.810 


1.3 


0.81 


27.36 


S98 


1111+1955 


0.299 


1.493 


17.7 


75.7 


<10.0 


<0.18 


>1360 


CI II 


G05 


1.152 





.648 


0. 


370 


0.4-1.4 


1.02 


26.79 


W92,G91,H07 


1120+1420 


0.362 


1.875 


85.06 


415.9 










S95 


2.460 


1 


.000 


0. 


613 


0.5 


0.8 


27.00 


S98 


1148+5254 


1.632 


12.228 


0.8 


4.3 










O06 


0.108 





.414 


0. 


450 


7.9 


-0.15(0.4) 


27.96 


O07 


1247+6723 


0.107 


0.479 


7.4 


13.9 


38.0 


0.26 


175 


cm 


P09 


0.344 





.191 


0. 


133 


0.4-1.4 


0.65 


25.06 


W92,G91,P92 


1256+5652 


0.042 


0.181 


70.0 


56.0 










VCS 


0.288 





.419 


0. 


.257 


<0.07 


0.89 


24.69 


W92,G91,H07 


1309+4047 


2.91 


24.905 


0.8 


6.3 








C1I 


OD12 


0.037 





.131 


0. 


115 


5.4 


0.23(0.8) 


28.22 


OD12 


1324+4048 


0.496 


2.719 


5.6 


32.7 


4.2 


0.12 


870 


CI II 


A12 


0.357 





.413 


0. 


246 


1.4-4.9 


0.94 


27.06 


W92,G91,P92 


1326+3154 


0.368 


1.910 


56.4 


278.1 


40.0 


0.88 


1030 


CI II 


K98 


4.631 


2 


.350 


1. 


630 


0.5 


0.66 


27.32 


S98 


1335+4542 


2.449 


20.183 


1.3 


4.1 










O06 







.821 


0. 


646 


5.0 


0.43(0.7) 


28.88 


O07 


1335+5844 


0.57 


3.31 


13.3 


84.2 


4.7 


0.16 


1800 


CI I 


A12 


0.299 





.744 


0. 


726 


6.0 


0.04(0.3) 


27.07 


O07 


1400+6210 


0.431 


2.300 


57.6 


313.3 










D95 


4.490 


1 


.800 


1. 


200 


0.5 


0.7 


27.41 


S98 


1407+2827 


0.0766 


0.333 


8.0 


11.1 


25.0 


0.13 


255 


CI I 


W12 


0.865 


2 


.532 


2. 


071 


5.1 


0.36(1.0) 


25.85 


O07 


1414+4554 


0.186 


0.874 


28.8 


86.7 


<14.0 


<0.16 


>1740 


CI I 


G05 


0.360 





.213 


0. 


139 


~1.4 


0.78 


25.69 


W92,G91,P92 


1415+1320 


0.247 


1.202 


8.1 


30.2 


60.0 


0.91 


100 


CD 


G05 


1.206 





.842 


1. 


564 


8.4-15.0 


-1.11(0.6) 


27.04 c 


W92,G91,H07 


1511+0518 


0.084 


0.370 


4.8 


7.3 


14.0 


0.08 


300 


CI I 


A12 


0.092 





.608 


0. 


811 


10.7 


-0.51(0.3) 


25.31 


O07 


1546+0026 


0.550 


3.087 


6.0 


37.4 










G05 


1.808 


1 


.149 


0. 


889 


<1.4 


0.47 


27.27 


W92,G91,H07 


1559+5924 


0.0602 


0.259 


7.5 


8.4 










B04 


0.229 





.197 


0. 


130 


0.4-1.4 


0.76 


24.59 


W92,G91,H07 


1609+2641 


0.473 


2.571 


50.6 


290.0 


23.0 


0.63 


1500 


CI I 


N06 


4.860 


1 


.710 


0. 


960 


1.0 


1.05 


27.71 


S98 


1616+0459 


3.197 


27.908 


1.4 


10.8 










O06 


0.317 





.771 


0. 


559 


5.0 


0.57(1.0) 


29.49 


D00,T05 


1723-6500 


0.0144 


0.060 


10.0 


2.8 


<106 


<0.10 


>90 


CI II 


PC03 


3.540 


4. 


.640 


3. 


670 


1.4-2.5 


0.40 


24.48 


T03 


1734+0926 


0.735 


4.402 


14.2 


100.5 


6.3 


0.25 


1300 


CI II 


A12 


1.110 





.740 


0. 


490 


2.3 


0.75 


27.58 


S98 


1755+6236 


0.027 


0.1148 


8.3 


4.4 










B04 


0.288 





.198 


0. 


147 


<0.3 


0.54 


23.78 


C98,G91,H07 


1816+3457 


0.245 


1.185 


36.0 


135.0 










G05 


0.678 





.355 


0. 


217 


~0.4 


0.90 


26.25 


W92,G05,G05 


1823+7938 


0.224 


1.072 


15.8 


54.8 










VCS 


0.297 






0. 


592 


-8.4 


(1.0) 


26.65 


W92,G05,G05 


1845+3541 


0.764 


4.615 


5.6 


31.0 


13.0 


0.55 


180 


CI I 


PC03 


1.031 





.794 


0. 


562 


1.4-4.9 


0.62 


27.55 


W92,G91,P92 


1939-6342 


0.181 


0.846 


43.4 


128.0 


26.0 


0.30 


1400 


CI II 


GP09 


14.98 


5 


.840 


2. 


425 


~1.4 


1.62 


27.63 


T03,H07 


1944+5448 


0.263 


1.285 


42.2 


165.0 


31.0 


0.50 


1080 


CI I 


PC03 


1.647 





.938 


0. 


610 


0.4-1.4 


0.78 


26.67 


W92,G91,P92 


1945+7055 


0.101 


0.445 


32.0 


57.0 


<30 


<0.19 


>1000 


CI I 


PC03 


0.921 





.645 


0. 


477 


0.4-1.4 


0.54 


25.46 


W92,G91,P92 


2022+6136 


0.227 


1.086 


7.0 


24.5 


13.0 


0.18 


440 


CI I 


PC03 


2.110 


2 


.820 


3. 


.210 


8.4 


-0.23(0.6) 


27.09 6 


S98 


2203+1007 


1.005 


6.488 


10.5 


81.8 


10.3 


0.53 


500 


CI I 


A12 


0.107 





.311 


0. 


231 


4.8 


0.53(1.0) 


27.86 


O07 


2355+4950 


0.238 


1.144 


49.3 


179.2 


21.1 


0.31 


1900 


CI I 


OCP99 


2.341 


1 


.552 


0. 


992 


0.7 


0.81 


26.81 


W92,G91,P92,S98 



*: Sources have both redshift and kinematic data and are used in Figures 3—10; 

a : Source structure identifications are: Class I and disturbed structure, Class II structure, and CD for core-dominated. One 
source has also been classified as a High-Frequency-Peaker; 

b : For CSOs with a spectral peak higher than 4.8 GHz, the spectral index is calculated from the optically thin part of the 
spectrum above the turnover frequency and given in bracket; 

c : The 1.4 GHz radio power is extrapolated from the 8.4-GHz flux density and the spectral index Og 5 ^/^- 
The references are : A 12: (lAn et al. 20T2l; B04: (IBondi et al. 200411; D95: jPalTacasa et al. 199511; F00: 
l IFomalont et al. 200011; G91 : I jGregory fc Condon 1991); G 03: flGiroIetti et al. 2003H; G05 I jGugliucci et al. 2"005); GP09: 
l|GiroIetti fc Polat idis 2009 ); K98: (|Kellermann et al. 1998); K04: (|KeUermann et al. 20"04); LOO: ijLiu et al. 2000) ; N06: 
I Nagai et al. 2006); O 06: llOrienti et al. 200611 ; O07: llOrienti et al. 2007): OC98: (jOwsianik fcTon way 1998); OCP99: 
jpwsianik et al. 19990 ; OD12 : jOrienti fc Dallacasa 2012jl : OM lO: l|Orienti et al. 20100 P09 : UPo latidis 20 09); P9 2: 
JPatnaik et al. 199211; PC 03: ^Polatidis fc Conway 20031 ; PT00: l|Peck fc T aylor 2000} ; S9 5: l[Sanghera et al. 1995) ; S98: 
l |StangheIIini et al. 1998); SOirptanghellini et al. 200l| ; T96: QTaylor et al. 19960; T03 : I jTingay & Kool 2003|>; VCS: VLBA 
Calibrator Survey UPetrov et al. 200811 . http://astrogeo.org/vcs/: W12: IIWu et al. 20121 : W92: dWhite fc Becker 1992L 
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Fig. 8. — Variation of spectral index with 
evolutionary time and with linear extent. 
The spectral index between 4.8 and 8.4 
GHz varies from a = 0.3 to unity for the 
range kinematic ages and of linear size. 
The four HFP sources with inverted spec- 
tra have a spectral turnover above 5 GHz. 
The source with the steepest spectrum is 
most likely a relic. 
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Fig. 9. — Histogram of kinematic ages for 
the sample of 24 CSOs with known red- 
shifts and expansion velocities from Ta- 
ble 1. The open squares represent lower 
limits to the kinematic age. 
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Fig. 10. — Histogram of the projected 
source size for a larger sample of 46 CSOs 
with known redshifts and size information 
from Table 1. 
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Fig. 1. — The radio power versus the linear extent of large scale radio sources: Compact, 
Medium-sized and Large Symmetric Objects. Exemplary evolutionary tracks based on para- 
metric modeling are depicted for the high-radio-power and low-radio-power sources using 
red and blue dashed lines. The black dashed lines marks the (approximate) boundary re- 
gion between stable laminar jet flows (above the lines) and unstable turbulent flows (below 
the lines). Symbols represent different morphological and spectral classes of radio sources 
(number of data points in bracket): black square: CSO (46), black circle: low-power GPS 
(12), red diamond: high-power GPS (27), purple cross: HFP (12), green circle: low-power 
CSS (114), blue open triangle: high-power CSS (44), blue filled triangle: FRII (140), and 
green filled star: FRI (79). Catalog references for Figure 1: Akujor & Garrington 1995 
(High-Power CSS); Black et al. 1992 (FR II); Bogers et al. 1994 (FR II); Bridle et al. 
1994 (FR II); Dallacasa et al. 2002a (High-Power CSS); Dallacasa et al. 2002b (High-Power 
CSS); de Vries et al. 2009 (Low-Power GPS); Fanti et al. 1987 (FR I and FR II); Fanti 
et al. 2001 (High-Power CSS); Kunert-Bajraszewska et al. 2010 (Low-Power CSS; K10); 
Kunert-Bajraszewska et al. 2006 (Low- Power CSS; K06); Laing et al. 1983 (FR I and FR 
II); Leahy et al. 1989 (FR II); Leahy & Perley 1991 (FR I and FR II); Liu et al. 2007 
(High-Power GPS; L07); Orienti & Dallacasa 2012 (HFP; OD12); Orienti et al. 2007 (HFP; 
O07); Orienti et al. 2010 (HFP; OD10); Orienti et al. 2006 (HFP; O06); Stanghellini et 
al. 1998 (High-Power GPS; S98) and Xiang et al. 2006 (High-Power GPS; L06). (A color 
version of this figure is available in the online journal.) 
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Fig. 2. — Evolutionary pathways for extragalactic radio sources. Mtype2 double sources 
will continue to grow in size and move horizontally across the diagram (panel a: red arrow). 
Low-power Mtype 1 double sources continuouly grow till a distance of about 100 kpc where 
the jet becomes unstable and flaring. Mtype 3 dying sources with decreasing or terminating 
power will move downward in the diagram and end up with relics (panel b). Re-energized 
and intermittent Mtype 4 re-started sources will move upward into the standard evolutionary 
tracks ( (panel c). Some low-power CSOs may evolve into Mtype 5 sources whose charac- 
teristic morphology is a hotpot within 1 kpc and a plume-like diffuse structure beyond the 
hotspot (panel d). Mtype 6 sources are shown as core-jet objects occupying the top-left 
corner in the P-D diagram. Some Mtype 1 sources, which are shaped by the ram pressure 
of the IGM, are classified as Mtype 7 (panel d). All sources below the jet stability line are 
not able to evolve into Mtype 2 double sources (panel a). 
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